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Abstract 

Using a (3+l)-dimensional solution of the relativistic Euler-equations for Pb + Pb at 160 
AGeV, space-time extensions of kaon emission zones are calculated from space-time densities and 
compared to the inverse widths of two- kaon Bose-Einstein correlation functions. The comparison 
shows a satisfactory agreement and it is concluded that because of the Gaussian shape of the 
kaon correlation functions, the space-time parameters of the kaon source can be calculated 
directly from space-time densities. In the case of intensity interferometry of identical pions this 
simplification is not recommended when applying Gaussian fits because of the present strong 
effects of resonance decays. The whole discussion is based on the assumption that hadron 
emission in ultra-relativistic heavy-ion collisions is purely chaotic or that coherence is at least 
negligible. 
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The knowledge about the lifetime and the spatial extensions of hadronic sources in ultrarelativistic 
nucleus-nucleus collisions constitutes further information in the search for a new state of nuclear 
matter, the quark-gluon plasma (QGP). In this effort an experimental technique has been developed 
which is dedicated to the measurement of Bose-Einstein correlations (BEC) jl]]. BEC functions of 
identical bosons are quantum-statistical observables, which in general contain information about 
ten quantities Q which characterize a hadron source: lifetimes, longitudinal and transverse exten- 
sions of the chaotic and the coherent source, temporal and spatial (longitudinal and transverse) 
coherence lengths and the chaoticty. In the case of pion interferometry it has been argued that a 
possible coherent subcomponent of the hadron source is almost inobservable || due to the appar- 
ently large contributions from resonance decays to BEC. Many parametrizations of experimental 
BEC functions are based on a neglect of a coherent component, resulting in the description of the 
hadron emitting source with only three quantities (e.g., lifetime, longitudinal and transverse radii). 
In ref. Q it was shown, that in fact even in the case of a purely chaotic hadron source three quan- 
tities are not adequate to fully parametrize the system under consideration. E.g., a cylindrically 
symmetric expanding fireball has in general to be described with an additional "cross-term" radius. 
Depending on the concepts of introducing a radius or a lifetime into a model of relativistically ex- 
panding hadron sources, it appears to be very likely that a direct relation between interferometry 
radii (better: inverse widths of BEC functions) and the true space-time characteristics of the source 
can be established. 

It is the purpose of this paper to study whether and how the true space-time features of hadron 
emitting sources are reflected in two-particle Bose-Einstein correlation functions of purely chaotic 
sources by applying a relativistic hydrodynamical description. Many hydro dynamical models (c/. 
H and refs. therein) are available which describe the dynamics of relativistic heavy-ion collisions. 
HYLANDER || belongs to the class of models which apply (3+l)-dimensional relativistic one- 
fluid-dynamics. It provides fully three-dimensional solutions of the hydrodynamical relativistic 
Euler-equations @. The model has been successfully applied to several different heavy ion colli- 
sions at SPS energies Q . Here some further results of the HYLANDER analysis [||, || of Pb + Pb 
at 160 AGeV (central collisions) will be presented. 
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From the hydro dynamical analysis [||, |J of Pb + Pb at 160 AGeV, it has been found that the 
K~ single inclusive spectra consist mainly of directly emitted kaons; less than 10% come from 
the K* resonance. Because the K* resonance contributions do not exhibit any significant effect in 
the inverse widths and the functional shapes of calculated K~K~ BEC functions, these resonance 
contributions have been neglected. Introducing the average and the relative four-momentum, K u = 
\{ki + A;^ ) and q u = k\ — k^, of two identical kaons emitted with their individual four- momenta, k\ 
and &2 j respectively, the Bose-Einstein correlation function of the kaon-pair is in terms of relativistic 
hydrodynamics defined as (cf. refs. ||, |l0|, O] ) 



Jd 4 x g d £{x u ,Kv)e iql,Xv j 2 
I d*x gf(x v , k\) ■ J d*x gf{x v , k%] 



where q° = E\ - E 2 , K° = (E 1 + E 2 )/2, and 

„dir(„ i„v\ _ 1 / k u d(T u (x' u ) 5 A (x u — x' u ) 
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In eq. (||) dcr^, and Ty are the differential volume element of the freeze-out hypersurface S, 
the 4- velocity of the fluid and the freeze-out temperature at space-time point x u , respectively. The 
quantity us is the strangeness chemical potential at space-time point x v . (For a calculation of BEC 
of identical pions including the decay of resonances, cf. refs. ||, [TT[| .) The function g^ r can be 
interpreted as the quantum analogue of the mean number of kaons of four-momentum k v at the 
space-time point x v , i.e., one can consider gj^ r as the quantum analogue of the space-time density 
distribution of the source emitting kaons of fixed four- momenta k v . For the following, let us keep 
the latter interpretation in mind. 



Fig. 1 shows a comparison of BEC functions of ir~ (solid lines) and of K~ (dashed lines) at K± = 
and K± = 1 GeV/c , respectively. The dotted lines correspond to the BEC functions of directly 
(thermally) produced tt~ . It can be seen that BEC functions of tt~ are strongly distorted by the 
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decay of resonances. In particular, the functional shape of the BEC functions of tt~ (solid lines) 
depends strongly on the momentum of the pair, i.e., it varies between the extremes of an exponen- 
tial and a Gaussian form. On the contrary, the functional shape of the BEC functions of thermally 
produced 7r~ and those of K~ is of almost perfect Gaussian shape, independent of the average 
momentum of the pair. In case of the shown transverse momentum-dependent BEC functions, the 
ones for thermal pions and kaons are almost identical. In the following we shall discuss only BEC 
functions of K~ pairs. 

Let us now require for a fit to a two-particle BEC function: 

(a) The selected functional shape of the fit should reproduce the entire two-particle BEC function. 

(b) The selected fitting function should respond to the changes in the average momentum of the 

pair. 



These demands ensure that a reconstruction of the original BEC function is possible through 
the knowledge of only the fitted parameters. In the case of kaon interferometry, (a) and (b) 
are easily to fulfill. The BEC functions can be fitted with a Gaussian, because their functional 
shape is independent of the choice of the average momentum of the kaon-pair under consideration. 
Therefore, we fit our results to (for the choice of the variables, cf., e.g. refs. |i~2]| ) 

C 2 (ki,k 2 ) = 1 + exp -q\Rj; ong {K) - q 2 ide R 2 ide (K) - q 2 out R 2 out {K) - 2 q\\ qout R 2 cross {K) . (3) 



A characteristic property of particle production from an expanding source is a correlation between 



the space-time point where a particle is emitted and its energy-momentum [13]. As a consequence, 
the inverse widths Ri(K) (i = ||, side, out, cross) extracted from Bose-Einstein correlation functions 
show a characteristic dependence of the average momentum of the pair, K {cf. also Fig. 1). 



Figs. 2a and 2b show the calculations for R S ide, Rout and Riong as functions of rapidity yx and 
transverse momentum K± from the fit of eq. (|3|) to the BEC functions in one (/-dimension, i.e., 



the remaining momentum differences have been set equal to zero. R cr oss has therefore not been 
extracted from a fit. For the specific choices of yx = and K± = 0, the inverse widths of the 
kaon BEC from Pb + Pb at 160 AGeV have already been shownQ in ref. |j. Since measurements 
of kaon interferometry radii have not been published yet for this particular heavy-ion reaction, the 
results shown in Figs. 2a and 2b have to be considered as further predictions of the hydrodynamical 
analysis presented in refs. [||, g|. 



In the following, we introduce an approximation [14, 15] of the original definition (c/. eq. ([I])) of 
the two-particle BEC, 



J d 4 x g$ r (x v ,K") e i q " x " \ 2 
Jd*x gf r {x v ,K» 



C2(h,h) * i + IJ ;;^r;; il, 2 1 - , (4) 



where q° = E\ — Ei and E% = K° = y w? + K 2 . In connection with analytical parametrizations of 
the momentum-dependent space-time density distribution gj^ r , the approximation eq. (Q) has been 
shown to be acceptable ||, |l5], |l6], 0. Here, we are going to apply the hydrodynamical solution 
already presented in refs. || |9| for g^ r according to eq. (|2|). 

Using x u = (t,x,y,z), we obtain for kaon pairs with \q\ <C Ek approximately 

q"x v « + - q ou tx - q si d e y ~ q\\z, (5) 



where = Ki/Ex {i = ||>-L). For any cylindrically symmetric emission function which can be 
roughly expressed in Gaussian form one finds with eq. @ when expanding exp[iq v x v ] in eq. @ 
for q u x v <C 1 the expression 

C 2 (K,q) « 1 + {1 - q 2 slde (y 2 ) ~ ([qout(P±t - x) + q^t - z)} 2 ) 

+ (qout((3±t - x) + q\\{(5\\t - z)) 2 + 0[{q u x v f)} , (6) 
1 Note the reduction of the inverse widths by a factor 1 / \pl. The inverse widths of ref. || have been obtained 
with a different parametrization. 
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where 

The average (£) can be considered as an expectation value from space-time densities g^ r - After 
exponentiating eq. (||), for any cylindrically symmetric system the BEC function C2 can be ex- 
pressed for small momentum differences {qiRi <C 1) in the form of eq. (]3|). When comparing the 
coefficients of the qiqj terms, the functions R 2 (K ) are found to be given through (c/. pi and refs. 
therein) 
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(3±t)(z - /3||t) = o- x - p ± t, z - P\\t ■ 


(8) 



Thus the inverse widths of the two-kaon BEC function can be calculated directly from the ten 
space-time averages (t), (x), (z), (t 2 ), {x 2 }, {y 2 }, (z 2 ), {xt), (zt) and (xz) (note that (y) = for 
a cylindrically symmetric expanding system). The quantities af and a^j in eqs. (||) can be in- 
terpreted as the variance (c/. also ref. [i~8[| ) of the random variable i and the covariance of the 
random variables i and j, respectively (of course, x v = (t, x, y, z) 6 E). The inverse widths of kaon 
BEC functions have therefore a geometrical interpretation involving the relativistic kinematics of 
the emitted particle. 



In Figs. 2a and 2b the space-time averages from space-time densities have been calculated from 
the r.h.s of eqs. (||) and are compared to the inverse widths of kaon BEC functions obtained from 
the direct fit (|3[). The agreement is remarkable, although R s ide an d Rout from fits to kaon BEC 
functions are always slightly underestimated, while Riong is always slightly overestimated. The 
differences of Ri on gi Rside an d Rout of both calculations have their origin in the approximations 
(H) " © an d in the fact that the kaon BEC functions are not of perfect Gaussian shape for all 
average momenta of the kaon pairs. The maximal differences are 0.54 fm, 0.15 fm and 0.17 fm, 
when comparing R[ ong , R s ide an d R ou t of both calculations, respectively. If one is willing to accept 
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errors of the magnitude here presented, one is tempted to avoid the calculation and subsequent 
fit of two-kaon BEC functions, since the inverse widths Ri(K) can be calculated in a much more 
effortless way directly from space-time averages. The feature of BEC of identical kaons to be of 
Gaussian shape leads through the knowledge of the functions Ri{K) always to an acceptable recon- 
struction of single BEC functions. In addition to the inverse widths which have been extracted from 
BEC functions of kaons, the cross-term radius R cross {K) has been calculated from the space-time 
expectation values (c/. Fig. 2b). 



Let us now return to the BEC functions of identical pions. We have found that the inverse widths 
of BEC of direct pions can be also treated as in the case of interferometry of identical kaons, i.e., 
through the use of effective Gaussians. The maximal differences when comparing Ri ong , R s ide an d 
Rout of both calculations are 0.78 fm, 0.20 fm and 0.20 fm, respectively. But it is obvious from 
the above presentation that an effective Gaussian would not constitute an appropriate fit for a 
single BEC function of all pions, i.e., including those from the decay of resonances, because of the 
drastic change in the functional shape (c/. Fig. 1). When it comes to the fit of experimentally ob- 
served BEC of pions, the situation is even more complicated, although measured BEC might have 
a Gaussian shape^J. Experimentally obtained BEC functions have to be calculated from separate 
integrations of the numerator and the denominator of the correlators of BEC functions of single 
momentum pairs with respect to detector acceptances ||, 19]. 



To summarize, Bose-Einstein correlation functions of identical kaons have been calculated for 
Pb + Pb at 160 AGeV from a (3+l)-dimensional solution of the relativistic Euler-equations. From 
the hydrodynamical treatment it has been found that almost all kaon BEC are of Gaussian shape, 
irrespective of the average momenta of the kaon pairs under consideration. Since in the calculations 
the kaon source was assumed to be completely chaotic, the BEC of kaons have been parametrized 
with a Gaussian and inverse widths of the BEC functions have been extracted. 
Interpreting the source function gj^ r (x I/ , k u ) as the quantum analogue to the space-time density 
2 It is therefore unclear, how the conclusions of ref. (l8j are affected. 
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distribution of emitted kaons of fixed four-momentum k v , space-time averages of space-time coor- 
dinates of the source have been calculated and are related to the inverse widths of BEC of kaons. 
The comparison of inverse widths of kaon BEC functions and space-time extensions from space- 
time densities shows a satisfactory agreement. It is concluded that because of the Gaussian shape 
of the kaon correlation functions, the space-time parameters of the kaon source can be calculated 
directly from space-time densities within errors of several 0.1 fm. 

In the case of intensity interferometry of identical pions this simplification is not recommended 
when applying Gaussian fits because of the presence of resonance decays the knowledge of only the 
inverse widths of the two-pion BEC is not sufficient to reconstruct the entire correlation function. 

B.R.S. is indebted to U. Heinz for discussions that initiated this work. Many helpful discussions 
with D. Strottman and his kind hospitality extended to the author at the Los Alamos National Lab- 
oratory are gratefully acknowledged. This work has been supported by the University of California 
and the Deutsche Forschunggemeinschaft (DFG). The author acknowledges a DFG postdoctoral 
fellowship. 
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Figure Captions 



Fig. 1 Examples of Bose-Einstein correlation functions of all n (solid lines), thermal it , and 
K~ (dashed lines), for K± = and K± = lGeV/c. The BEC functions have been calculated 
for Pb + Pb at 160 AGeV. 

Fig. 2a Momentum-dependence of R s id e and R ou t extracted from Bose-Einstein correlation func- 
tions of kaons compared to the directly from space-time densities calculated ones. 

Fig. 2b Momentum-dependence of Ri ong extracted from Bose-Einstein correlation functions of 
kaons (BEC) compared to the directly from space-time densities calculated ones. The radius 
Rcross has been calculated from space-time densities only. 
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Figure 1 
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Figure 2a 
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Figure 2b 



